Large spin relaxation rates in trapped submerged-shell atoms 
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Abstract 

Spin relaxation due to atom-atom collisions is measured for magnetically trapped erbium and 
thulium atoms at a temperature near 500 mK. The rate constants for Er-Er and Tm-Tm col- 
lisions are 3.0 x 10^^'^ cm^ s~^ and 1.1 x 10^^'^ cm'^ s^^, respectively, 2-3 orders of magnitude 
larger than those observed for highly magnetic 5-state atoms. This is strong evidence for an addi- 
tional, dominant, spin relaxation mechanism, electrostatic anisotropy, in collisions between these 
"submerged-shell" L 7^ atoms. These large spin relaxation rates imply that evaporative cooling 
of these atoms in a magnetic trap will be highly inefficient. 

PACS numbers: 34.50.-s,37.10.De,32.70.Jz 
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Research in cold and ultracold atoms has in recent years increasingly broadened scope 
beyond the alkali atoms to explore and exploit the diverse range of atomic and chemical 
properties found across the periodic table. In particular, the lanthanide rare-earth (RE) 
atoms have attracted considerable experimental and theoretical interest. Recent experiments 
with RE atoms have resulted in Bose-Einstein condensation of Yb magneto-optical 
trapping of Er \^ and Dy [s], Zeeman slowing of Tm {3], and large ensembles (> 10^^ 
atorns) of buffer-gas loaded and magnetically trapped RE atoms of several species below 
1 K [sl, l6|. This interest in RE systems stems from important, sometimes unique, attributes 
such as narrow transitions which allow for low Doppler cooling limits and improved frequency 
standards 0,0], large magnetic moments with strong long-range dipolar interactions, and 
"submerged-shell" character that in certain circumstances can shield atom-atom interactions 
from anisotropic valence electron shells . Progress with these systems — or any novel atomic 
system — is dependent on coUisional processes, in particular low rates of inelastic collisions 
including spin relaxation collisions in trapped samples. Spin relaxation can cause heating 
as well as drive atoms out of the desired quantum state, thus preventing cooling to lower 
temperatures and limiting experimental sensitivity and the capacity for new discovery. 

Previous experiments and theoretical work with RE atoms, including Er and Tm, dis- 
covered suppression of electrostatic anisotropy in RE-helium collisions. Specifically, in this 
interaction the anisotropic 4f electron distribution was found to be shielded by closed 5s 
and 6s electron shells This "submerged-shell" nature allowed for sympa- 

thetic cooling of RE atoms by cold He and efficient buffer-gas trapping of large numbers of 
atoms (> 10^^) at millikelvin temperatures. It also explained reduced coUisional frequency 
broadening of hyperfine clock transitions in Tm [13j. The discovery of efficient shielding 
and consequent low inelastic rates in the RE-He system gave hope that similar suppression 
would be found in RE-RE collisions and could allow for efficient evaporative cooling fl^ . 
This could, for example, provide a path to quantum degeneracy for magnetically trapped 
RE atoms. 

In this Letter, we present measurements of spin relaxation rates in two-body collisions 
of trapped submerged-shell species Er {L = 5) and Tm (L = 3), finding them to be very 
large, in striking contrast to the low rates observed in RE-He systems. These large rates 
imply an additional spin relaxation mechanism other than spin exchange, dipolar relaxation, 
and second-order spin-orbit coupling, which are well-known from studies with alkali atoms. 
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Although theory has proven very effective for understanding the RE-He system, the current 
theory of RE-RE cold collisions is incomplete. Despite the development by Krems et al. of a 
theoretical framework for collisions of two L 7^ atoms 15| , in the RE case these calculations 
are extremely difficult and to our knowledge no theoretical predictions yet exist. Recently, 



an experiment was done with the transition metal titanium at a temperature of 5 K 16|. 

n n 

Ti has submerged-shell structure, as observed in Ti-He collisions [17|, 118|, but an order of 



magnitude weaker shielding of electrostatic anisotropy [10|. Rapid decay of ^°Ti electron 
spin polarization was observed due to ^°Ti-Ti collisions, but the mechanism of this loss 
could not be determined because spin relaxation could not be separated from spin exchange 
with unpolarized Ti isotopes. Thus, whether or not submerged-shell atoms exhibit low spin 
relaxation rates — indicating submerged-shell behavior in these processes — had remained an 
open question. 

Our experiment takes place in a double-walled plastic cell maintained at a temperature 
of ~ 500 mK by a superfiuid helium heat link to a dilution refrigerator (see Fig. [1]). We 
produce either trapped atomic Er or Tm by laser ablation of solid metal foils into ^He 
buffer-gas in the presence of a magnetic quadrupole field (trap depth up to 3.7 T) produced 
by large superconducting anti-Helmholtz coils surrounding the cell. The ablated atoms cool 
via elastic collisions with the cold buffer-gas and within 50 ms [19] assume a Boltzmann 
distribution in the trap with peak density of up to 7 x 10^^ cm~^. The trapped cloud is 
interrogated via laser absorption spectroscopy on the 400.9 nm (J = 6 ^ 7) and 415.2 nm 
( J = 6 ^ 5) transitions of Er and on the 409.5 nm (J = 7/2 5/2) transition of Tm. 

The amount of buffer-gas in the cell is regulated such that the He density is sufficient 
to cool the atoms after ablation, but insufficient to cause significant atom loss from RE-He 
collisions. This regulation is accomplished by independent control of the cell temperature 
and the amount of He initially present in the cell. The lack of observed loss from buffer-gas 
collisions 1 second after ablation implies a He density less than 10^^ cm~'^ [191]. Since the 
observed ablation yield implies a higher initial buffer-gas density, it is likely that heating 
from the ~ 5 mJ ablation pulse temporarily desorbs additional He from the cell walls, which 
re-adsorbs rapidly. We deliberately maintain a buffer-gas density of ~ 10^^ cm~'^ after trap 
loading in order to maintain thermal equilibrium between the trapped atoms and the cell. 

Example spectra of magnetically trapped Er and Tm are shown in Fig. [21 showing peaks 
for both Zeeman-broadened Amj = ±1 and narrow Arrij = transitions. The relative 
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magnitudes of spectral features may be used to estimate the mj state distribution, however 
for the case of a probe laser passing through the trap center the absorption is primarily 
determined by the total peak atom density rather than the contributions from individual 
states. Isotope shifts for the 400.9 nm transition of Er were not found in the literature and 
were determined for nuclear spin-zero isotopes by fitting to spectra measured in zero field at 
~ 4 K. The shifts for isotopes ^^^Er, ^^^^Er, and ^^°Er from the ^^*^Er peak are -0.80(4) GHz, 
0.81(1) GHz, and 1.66(2) GHz, respectively. 

As noted above, we ensure that the He density is sufficiently low such that neither elastic 
nor inelastic RE-He coUisional loss is observed (see Fig. [3]) . The trap loss is then determined 
by the rate equation: 

n(f, t) = -[fevap{Etrap, T)gel + gin]n{r, tf , (l) 

where n is the local density of trapped atoms, and gei and gin are the rate constants for elastic 
and inelastic atom-atom collisions. The function fevap is the fraction of elastic collisions at 
temperature T that are energetic enough to produce atoms with energy above the trap depth 
Etrap such that the atoms will adsorb on the cold cell walls and be lost from the trap. In 
our experiments T is low enough such that fevap < 1% |l4) ^'^d thus elastic collisions do not 
significantly contribute to atom loss. Ignoring the first term in Eqn. [1], we solve for n{f,t), 
spatially integrate over the trap distribution, and take the reciprocal to reach the simple 
two-body decay result: 

+ (2) 



no{t) n(r = 0,t) no{t = 0) 

Plotting versus time yields a straight line of slope gin/8. Data for Er decay is plotted in 
this manner in Fig. [3] and fit to Eqn. [2l Additionally, a combined fit with free parameters for 
Er-Er and Er-He coUisional loss processes yields a Er-He decay rate consistent with zero, 
and therefore we conclude that the loss is indeed due to Er-Er collisions. 

Fits of atom loss to Eqn. [2] yield gin to be 1.5 ±0.2 x 10"^° cm^ s"^ for Er and 5.7 ± 1.5 x 
10^^^ cm^ s~^ for Tm, with accuracy limited by the density calibration determined from 
spectra. Both rates are significantly higher than inelastic rates observed for highly magnetic 
S-state atoms such as Cr, Eu, Mn, and Mo 

m 
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231]. The spin relaxation rate 



constants for these species were measured in similar magnetic traps at similar temperatures 
and found to be < 10"^^ cm'^ s~\ consistent with the magnetic dipole-dipole interaction 
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25| described by 



V,i,oie{^ = ^^[(Ji- J2) -3(Jrf)(J2-f)]. (3) 

47r 

where /i is the magnetic moment and r is the distance between two atoms with angular 
momenta Ji and J2, respectively. The spin relaxation rate constant is dependent on the 
specific form of the interatomic potential, however the general /i^ scaling implied by Eqn. [3] 
provides a relation between dipole-induced inelastic loss rates for atoms of similar electronic 
structure. Eu (/i = 7 fis) and Mn (/i = 5 fis) in particlar have submerged-shell character 
similar to Er (/x = Jub) and Tm (/x = 4 hb) [26|. Scaling the cross sections measured 
for Eu ^ and Mn j23| by /x"^ and averaging yields = 3.4 x 10~^^ cm^ s~^ for Er and 
Qin = 3.5 X 10"^*^ cm^ s~^ for Tm. The observed inelastic rate constants for Er and Tm 
in our experiments are 2-3 orders of magnitude larger than these scaled dipolar values, 
inconsistent with the dipolar loss model and implying another loss mechanism. 

A significant fraction of atoms (> 20%) in our experiments have mj ^ J, as determined 
by observing Amj = transitions on the 415.2 nm (J = 6 ^ 5) line of Er. Two-body 
electronic spin exchange collisions will tend to purify the atomic ensemble towards the 
mj = J state, but the stability of spectral features with time implies that this is not the 
case. In addition, since such collisions conserve the total mj, they cannot cause loss to 
untrapped states without also populating more strongly-trapped states, which would cause 
an unobserved net increase in absorption. Nuclear spin exchange could lead to trap loss, 
but observed rates for this process in other submerged-shell atoms with only / > isotopes 
have shown it to be much slower than the loss observed here [s], [3] ■ Hence the observed 
loss is spin relaxation to untrapped states. In our analysis, we assume gin to be the same 
for all pairs of atoms of any vtij. 

For spin relaxation collisions resulting in a final state with mj > 0, relaxation may not 
lead immediately to trap loss. In that case, the deduced from loss may be smaller 
than the true spin relaxation collision rate constant, which we will call Qst- Calculations for 
collisions between He and L 7^ atoms such as Tm and O yield larger rates for Amj = ±1,2 
transitions than for other transitions, creating effective selection rules jlll . 27 1. Although 



RE-RE sytems are not theoretically well-understood, if such selection rules held in the case 
of Er, the mj = J = 6 state would on average require several inelastic collisions to reach 
an untrapped state, contributing to the nonzero rrij ^ J state population noted above. In 
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addition, collisional energy can promote inelastically colliding atoms to higher mj states 
and inhibit loss. These thermal excitations are suppressed for gjfisB ^ kT, however this 
condition fails near the trap center where B = 0. Considering both these effects, the observed 
stability of the spectrum suggests that the rrij state distribution achieves a slowly-varying 
balance between loss and excitation. We confirmed this model with simulations of inelastic 
decay, including thermal excitations and exploring a range of initial mj state distributions 
and selection rules. The simulations suggest a ratio gsr/din of 2.0 tl'^- 

Currently there exist no theoretical predictions for RE-RE spin relaxation rates due to 
the complexity of RE electronic structure; however, one reasonable hypothesis to explain 
rapid spin relaxation of Er and Tm is that it is induced by electrostatic anisotropy, as 
is observed in anisotropic outer-shell systems. Experiments with metastable states of 
Ca and Yb have measured inelastic collision rate constants greater than 10~^^ cm'^ s~^ 



28|, |29|, nearly as large as the Ca*-Ca* and Yb*-Yb* elastic rate constants. These inelastic 



rates are similar to those we observe here for Er and Tm atom-atom collisions, suggestive 
of a complete lack of suppression of electrostatic anisotropy-driven spin relaxation and in 
contrast to the dramatic suppression of > 10^ observed for collisions with He. In the RE- 
He case, the Born-Oppenheimer potentials corresponding to different projections of the 
electronic angular momentum onto the internuclear axis are nearly degenerate [lOj, which 
leads to a low probability for reorientation of the magnetic moment. The RE-RE interaction 



potentials are much deeper than those of the RE-He system |30|, l3l| , and so it is possible 
that electrostatic anisotropy has a similarly stronger effect. 

In conclusion, we have measured the loss rate constants for inelastic Er-Er and Tm- 
Tm collisions and found them to be large. For comparison, the maximum elastic cross 
section agi in the absence of shape resonances can be derived from the well-known unitarity 
limit 



32|. Using the Ce constant calculated for the Yb-Yb system 



30( 1 and assuming elastic 



collisions between submerged-shell lanthanide RE atoms to be similar, we find the maximum 
del = <^eiv ~ 8 X lO^^'' cm^ s^^ at 500 mK. Hence the ratio Qei/ Qsr ^ 10 for both Er and Tm, 
implying that evaporative cooling of these atoms in a magnetic trap will be highly inefficient 



14| . At this temperature we expect ~ 40 partial waves to contribute to collisions, and we 



note that may be different in the ultracold s-wave limit. However, this limit is rather 
low for these heavy colliding atoms (10-100 /iK), so the multi-partial wave physics will be 
applicable over a range of experimental conditions. 
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The large spin relaxation rates for Er and Tm reported here, along with those reported 



33| . represent significant evidence that 



for Ti and recently measured separately for Dy 
the submerged-shell character exhibited by roughly a third of the periodic table and which is 
responsible for dramatic suppression effects in atom-He collisions does not imply suppression 
of electrostatic anisotropy-driven spin relaxation in atom-atom collisions. As a result, the 
highly successful method of evaporative cooling in a magnetic trap may remain confined to 
(isotropic) S-state atoms. In addition, lifetimes for optically trapped atoms in L 7^ states 
may be short due to spin relaxation unless trapped in the absolute ground state. 

We would like to acknowledge helpful discussions with Timur Tscherbul. This work was 
supported by the NSF under grant number 0757157 and through the Harvard/MIT Center 
for Ultracold Atoms. 
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FIG. 1: (color online) Diagram of the buffer-gas trapping apparatus. A cryogenic valve separates 
the trapping region of the experimental cell from an additional pumping region and can be used 
to regulate the buffer-gas density. 
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FIG. 2: (color online) (a) Absorption spectrum of Er on the 400.9 nm (J = 6 ^ 7) transition 
in a 0.99 T (4.6 K) deep magnetic trap at 530 mK with peak density of 4.6 x lO^'' cm~^. The 
Am J = +1 magnetically broadened peaks of the dominant isotopes are labeled. The sharper peaks 
are Amj = transitions. Hyperfine constants are unknown for the ^^"^Er isotope (23% abundance), 
and it is ignored in the spectrum simulation. Due to the substantial Zeeman broadening, this does 
not significantly affect the implied atom density and temperature, (b) Absorption spectrum of Tm 
on the 409.5 nm (J = 7/2 5/2) transition in a 3.3 T (8.8 K) deep trap at 500 mK with peak 
density of 3.8 x 10^^ cm~^. Tm has a single isotope with 1=1/2 and known hyperfine splitting 
34|. 
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FIG. 3: (color online) Er decay at 530 mK in a 0.99 T (4.6 K) deep magnetic trap after ablation 
at t = s. The vertical axis is the reciprocal of the peak atom density obtained from spectra. The 
solid red line is a fit to Eqn. [2l The dashed blue curve is a fit to the exponential decay expected 
for collisions with a constant He background. The excellent fit to Eqn.[2](r = 0.998) indicates that 
the atom loss is from Er-Er collisions. 
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